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A perturbation approach based on the first-order mean-spherical approximation (FMSA) is proposed. It con- 
sists in adopting a hard-sphere plus short-range attractive Yukawa fluid as the novel reference system, over 
which the perturbative solution of the Ornstein-Zernike equation is performed. A choice of the optimal range 
of the reference attraction is discussed. The results are compared against conventional FMSA/HS theory and 
Monte-Carlo simulation data for compressibility factor and vapor-liquid phase diagrams of the medium-ranged 
Yukawa fluid. The proposed theory keeps the same level of simplicity and transparency as the conventional 
FMSA/HS approach does, but turns out to be more accurate. 
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1. Introduction 

The most common and universal method in theories of fluids is the perturbation approach. 
Nonetheless, without analytic results of the integral equation theory for simple systems the per- 
turbation expansion could not be accomplished. To be specific, the analytic results of the mean- 
spherical approximation (MSA) for the fluid of hard spheres (for which the MSA is identical to the 
Percus-Yevick theory) [TJ [5] and the Yukawa fluid |3l H] are used to describe the properties of the 
reference system. In addition to these conventional MSA results there is also a perturbed version 
of the MSA called the first-order mean-spherical approximation (FMSA) [5] . 

When compared to computer simulation data for the Yukawa fluid model, the FMSA theory 
has been shown to be reasonably good but still less accurate than the common MSA [6J. However, 
an advantage of the FMSA theory, hereinafter referred to as FMSA/HS since it employs the hard 
sphere (HS) fluid for the reference, is its simplicity and transparency. These are the factors causing 
its increasing attractiveness and potential applicability in the liquid state theory and also, in 
general, in the soft condensed matter theory. The most important ideas of the FMSA/HS theory 
are outlined in the following section. 

The goal of this paper is to present a modification of the FMSA/HS theory. The modification 
pursues the idea that there is an alternative to the HS reference, namely a mean field theory based 
on a Yukawa fluid reference [7]. In section 3 the possible choices in this direction are considered. 
The results and discussions are the subject of section 4. The study is summarized in section 5. 

2. Conventional FMSA/HS theory 

The FMSA/HS is a theory based on the solution of the Ornstein-Zernike (OZ) equation. 



h{k) = c{k)+ ph{k)d{k) 



(1) 
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by employing the perturbation expansions limited to the first order correction term, 

h{k) = ho{k) + Ah{k) , 
c{k) = 5o(fc) + A5(fc), 

Q(fc) = Qo W + AQ (fc) , (2) 

for the total, h, and direct, c, correlation functions combined with Baxter factorization function, Q, 
respectively. Here and in what follows all symbols with the tilde will denote the three-dimensional 
Fourier transforms, while all symbols with the hats will denote the one-dimensional Fourier trans- 
forms or the Laplace transforms. The terms with subscript "o" denote the contribution of a reference 
system, and A/i, Ac, AQ are the respective corrections. 

The FMSA/HS theory has been developed by Tang and Lu [S HI [9] for the hard-core based 
fluid models 

/ X / oo, r < cr, , , 

[ 0(r), r > (T, 

with a being the diameter of a hard core and 4'{r) being the potential function for out-of-core pair 
interaction between molecules. In their studies Tang and Lu were considering the fluid of hard 
spheres (HS) as the reference (or unperturbed) system for which all three functions that enter 
expansions ([2]) are well known. In particular, for the total correlation function of the HS reference 
system, ho = /irs , Tang and Lu used the expression 

9o (s) = 5HS (s) = .2 A , , , > (4) 

that is the Laplace transform of the radial distribution function .gHs('') = ^hs(^) + 1 resulting from 
the solution of the OZ equation for HS fluid in the Percus-Yevick approximation [1]. The function 
Qhs (t) refers to the Laplace transform of the Baxter factorization function of the HS fluid, 

Qo [t) = Qm [tj = -2-: , (5j 

(1-77) t'i 

where 

S{t) = {l~Tjft^ + 6'n{l-Tj)f + 18r]h-12Tj{l + 2'n), (6) 
L{t) = (l + |)i+l + 2r; (7) 

and 77 = irpa^ /6 is the packing fraction. 

The main result obtained by Tang and Lu concerns an expression for correction Ah to the total 
correlation function (or radial distribution function) resulting from the out-of-core attraction (t>{r). 
It reads [5J 

P {ika) 
where 

00 

2Q|s(-ifca) 2m J {y-k)Ql^[-iya) 

— 00 

with function /SJJ (k) defined as 

00 

AU (fc) = J rAc (r) e-''='^dr. (10) 
In accordance with this approach the necessary closure to be used in equation (10) reads as 



Ahik) = ,^-^, (8) 



^^FMSA/HS ^ ^ fQj. ^ ^ 
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with j3 = l/k^T and T being the temperature. 

In general, without any relation to the above version of the FMSA theory, a key point of 
its application to particular fluid model [determined by the out-of-core interaction potential (t>{r)\ 
concerns an evaluation of the integral in equation (|9|. So far, within the framework of the FMSA/HS 
theory this has been done for the Yukawa, Lennard- Jones, Kihara and sticky fluids by Tang and 
Lu [5] and for the square- well fluids by Tang and Lu [9 , and by Hlushak et al. [TU] . In the particular 
case of the Yukawa (Y) fluid model, 

-z(r-cr) 

(r) = -ea , (12) 

r 

for which the function AU (k) ~ e"'''"", an integration contour in the right-hand side of equation ^ 
can be closed in the upper complex half-plane and evaluation of the integral is rather simple and 
does not require the calculation of the residues at zeroes of the function Qhs (— ifccr)- Then, the 
Laplace transform of the correction term ^ for the total correlation function of the Yukawa fluid 
within the FMSA/HS theory reads 

A/i™SA/HS (^) ^ . (13) 

(s -I- z) Q|s (s^) Qhs i^^) 

Once the radial distribution function g — go + Ah is known, the thermodynamics of the system 
can be calculated. 



3. FMSA theory based on Yukawa reference system 

The first order mean-spherical approximation theory that we are dealing with is a kind of 
perturbation theory approach. Usually, within the perturbation theory, in order to improve the 
performance of the first order approximation one should calculate the second order correction term. 
As an alternative to this common way here we propose to qualitatively modify the reference system 
over which the perturbation is calculated, and as a result to continue working within the same first 
order approximation. This is of particular importance since our aim is to keep the simplicity and 
transparency of the improved FMSA theory at the level established by the conventional FMSA /HS 
approach. 



3.1. Division of the potential 

Similar to the classical perturbation theory approach we proceed by dividing the initial inter- 
action potential u{r) into two parts 

u (r) = Uo (r) + Alt (r) , (14) 

that are the reference and residual contributions, respectively. However, for the reference system 
potential Uo (r) we will require that besides the hard-core repulsion it should (i) include a piece of 
the attractive tail of the same strength e as the total potential, and (ii) extend over the range that 
is somewhat shorter than the range of the total interaction. The function that allows us to satisfy 
and control these requirements in an easy and a natural way is the Yukawa (YO) potential. Thus, 
we may define the desired reference fiuid as follows: 

Mo (r) = UYo{r) ^ \ _ , r <a; ^^^^ 

I — ecre °^ > /r, r ^ a . 

The improved FMSA closure, referred to as FMSA/YO, reads as follows: 

^^FMSA/Yo ^ _ ;3eae-^°('-'^) /r, for r > tj , (16) 

where is the so far undefined parameter. 
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3.2. Yukawa reference system 

We note here that some authors have already suggested to utiHze the non-HS reference system 
in Hquid state theory (e.g., see |1H IT^). Moreover, the short-range attractive Yukawa fluid has 
been already considered as an alternative to the HS reference system by Melnyk et al. [7l [T3l - fl6] 
in their studies within the framework of an augmented van der Waals theory for simple fluids. 

It is very important that, like in the case of the HS reference potential, all properties for the 
Yukawa reference potential ( |15[ ) , including the Baxter function Qa = Qyo , are available in the 
literature. They can be obtained either within the MSA theory for Yukawa fluid [e.g., see Blum 
and Hoye Kalyuzhnyi et al. |17| fT5]] or within the conventional FMSA/HS theory due to Tang 
and Lu as it is described in previous section 2. The MSA theory is more accurate than FMSA/HS 
theory, but the latter is simpler. For the purpose of the present study we decided to sacrifice the 
accuracy in order to maintain simplicity and transparency. Thus, in what follows the FMSA/HS 
theory will be employed to describe properties of the YO reference system. And, as we will see 
hereinafter, despite a less accurate description of the reference system, the FMSA/YO theory still 
shows a notable improvement against the FMSA/HS theory. 

Following the conventional FMSA/HS approach, the Laplace transform of the radial distribution 
function of the YO reference system reads 

^ecre"'*'^ 

9o (s) = ffYo (s; 2o) = 5HS (s) + , — , , r , (17) 

Whs v^oCj 

where the first term corresponds to the contribution of a hard-sphere repulsion, while the second 
one is the contribution due to the short-range attraction attributed to the YO reference system in 



accordance with equation (151. Similarly, the FMSA/HS result for the Baxter factorization function 
of the YO reference is 

ZoCrQus {-t) c"* - {t + Zocr) Qm (*) 



Qo it) = Qyo (t; zoa) = Qhs (t) + 12/3£y? ^ / ' ° ^^""^ ^ . (18) 

tZa<7 [t + ZqCt) Q^g [ZoCr) 

3.3. Full system within the FMSA/YO theory 



After all properties of interest for the YO reference fiuid are specified, we turn to the entire 
system determined by the interaction potential u{r) or more precisely by the out-of-core potential 
function (/)(r). Although this function can be substituted by any potential function used in literature 
to represent the simple fluids [e.g., Lennard- Jones, Sutherland, Yukawa or Kihara potentials, etc.[ 
for the purpose of present study we proceed with the Yukawa (Y) fluid model already deflned 



according to equation (12 1 



By introducing Yukawa potential (pyir) into the FMSA/YO closure (16 1 we note that for both 
functions that enter Ac^^^*/^°(r), the function AU (k) is proportional to e~''"'°' and the integral 
in equation ([£]) can be easily evaluated in the way it was already discussed in Introduction section. 
Then, the correction term ^h^^SA/Yo ^i^^i jg necessary to evaluate the radial distribution function 
of the full system, 

g (s; Zo, z) = go (s; Zq) + A/i (s; Zo, z) , (19) 

is given by 



;.s + Zo] Wyo (.'Sc^; ^o'^) Wyo i^o^; ^o^) 

(20) 



(.s + Zo) Qyo ('Sc^; ^ofy) Qyo (^q^! ^o<y) 
(s -I- z) Q^r^ (scr; Zo(j) Q'^Q {za; ZoCr) 



The calculations of thermodynamics within the FMSA/YO approach can be made through the 
energy route, i.e., in the way that is quite similar to how it was done by Tang et al. within the 
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conventional FMSA/HS theory [6J. It consists in evaluating the internal energy using its definition 
through the radial distribution function, 

CO 

U f 

2npf3 / drr'^g (r; z^, z) u (r) = 12r]l3eae"^g {z; Zo, z) . (21) 



NkT 



This result is used to derive the Helmholtz free energy in the form 

A - Aid 



NkT ao + oi + 02 , (22) 



where 



4?7 - 3??' , . 

ao = -; ^2", (23) 

(1 - r/) 

BeL (za) ,„ 

(1 - 1]) QhS (zcr) (zcr) 



02 = -Gry/J^e^ 



(zcr + ZoCr) Q^g (zcr) Q^g (zqO-) 



(zcr + Zocr) Qyo (^'^^ ^o^^) Qyo (^o'^^ ^o^) ^zaQyo {za; Zocr) Qyo (^'^'^ ^o^) 



(25) 



In the limit Zq — z the above formulas reduce to the conventional FMSA/HS theory by Tang 
and Lu [6 . The chemical potential and pressure are obtained employing standard thermodynamic 
relations. 



4. Results and discussions 



4.1. How much of the attraction should be treated as the reference? 

The above outlined FMSA/YO theory con- 
tains the parameter Zo that determines the 
range of reference attraction and needs to be 
specified. The most straightforward approach 
to proceed deals with the minimization of the 



free energy ( 22 1 of the Yukawa system over the 



range of possible values of Zo . Noting that z^ 
enters only 02 contribution, we deduce that the 
outcome of the minimization, i.e., Zo,min , de- 
pends only on the density of the fluid but not 
on the temperature. The corresponding density 
dependence of Zo.min for several Yukawa fluids 
is presented in figure [T] It is evident that in the 
limit of a vanishing density and for all systems 
considered, the parameter Zo,min tends to the 
value of 2z. 

This rather remarkable result for choosing 
the YO reference system in such a strict form 
is quite probably limited to the FMSA theory 
only. Nevertheless, it is quite consistent with 
our earlier findings for the YO reference decay 
parameter within the framework of the Yukawa 
based van der Waals theory for the simple fluids 




Figure 1. Density dependence of the YO refer- 
ence decay parameter Zo.min that minimizes the 
free energy for Yukawa fluid models with differ- 
ent value of the decay exponent z as it is spec- 
ified for each curve. One can see that at low 
density Zo.min — >■ 2z. 

[71 [TBl - fTB] . First of all, there is a requirement for 
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-2^0, mill to be larger than z providing in this way the range of the reference attraction which will be 
shorter than that in the parent fluid. Secondly, it is obvious that the reference fluid should be in 
a one-phase region for the set of density and temperature parameters used in the studies of the 
parent fluid. In the case of a hard sphere reference, this requirement is satisfled since HS fluid does 
not exhibit the liquid/vapor transition. By adding an attraction one provides the possibility for 
the phase coexistence. However, the critical temperature is always getting lower if the range of 
attraction is shorter. In [14 we presented a collection of the Monte Carlo generated liquid/vapor 
envelopes for Yukawa fluid with different values of the decay parameter (see figure 1 in [14J) from 
which it follows that the values 3 < ZqU < 6 can be used as a reference decay for the Lennard- 
Jones-like Yukawa fiuid (za = 1.8) since the critical point temperature in such a reference fluid 
will always be lower than the triple point temperature in parent fluid. 

On the other hand, analyzing the dependence of the critical point coordinates of the Lennard- 
Jones-like Yukawa fluid (za = 1.8) on the reference system decay parameter Zq , in |13j we have 
shown that the best agreement with computer simulation data can be reached if 3 < ZqCt < 4. 
Thus, the result Zq — 2z agrees with these findings and in all our subsequent calculations for the 
Lennard- Jones-like Yukawa fluid we impose ZqCt = 3.6, making it independent of the fluid density. 

4.2. Comparison between FMSA/HS and FMSA/YO solutions 

To illustrate the improvements that the FMSA/YO theory yields over the conventional 
FMSA/HS theory, we compare the predictions made up by these two theoretical approaches for 
compressibility factor and vapor/liquid phase diagram of the most popular Yukawa fluid model 
defined by the decay parameter za = 1.8. 




Figure 2. Isotherms T* = 1.5 and T* = 1 

(from the top to the bottom) for compressibil- 
ity factor, Z = PV/NkT, of the Yukawa fluid 
with za = 1.8 as obtained from the conven- 
tional FMSA/HS theory (black dashed lines) 
and from the improved FMSA/YO theory with 
Zocr — 3.6 (red solid lines). The symbols denote 
computer simulation data by Schukia |19| . 



Figure 3. Vapor-liquid phase diagram of the 
Yukawa fluid with za = 1.8. The notations 
are the same as in figure [2] 



The compressibility factors, Z — PV/NkT, of the Yukawa fiuid with za = 1.8 and for two 
temperatures, T* — 1.5 and 1, that result from both FMSA approaches are presented in figure |2] 
To estimate the effectiveness of the theory, the computer simulation data due to Schukia [T^] are 
shown as well. As it was already pointed out by Tang [6], the conventional FMSA/HS theory does 
a good job in predicting the compressibility factor of the Yukawa fluid. Nevertheless, the proposed 
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FMSA/YO theory shows to be even more accurate, especially at lower temperature and in the 
region of intermediate densities. 

While the enhancement may seem to be minor in the case of compressibility factor, it is more 
pronounced for the vapor/luqiud phase diagram shown in figure [I The results of the FMSA/YO 
theory lead to the shrinkage of coexisting densities envelope in the region near the critical point, 
approaching the computer simulation data. The estimated improvement in comparison with the 
FMSA/HS predictions is around 30%. 



5. Conclusions 



Being compared with computer simulation data for the Lennard- Jones-like Yukawa fluid (za = 
1.8), the first order mean-spherical approximation (FMSA/HS) due to Tang et Lu has shown to 
be only slightly less accurate in the calculations of the thermodynamics and liquid/ vapor phase 
coexistence than the full MSA theory (details of this comparison and a corresponding discussion can 
be found in the |3]). At the same time, the FMSA/HS is a much simpler theory. Abbreviation HS 
underlines here that the theory is based on the hard-sphere reference system. In the present study 
we reported a further improvement of this approach, which we called a conventional FMSA/HS 
theory by introducing the short-ranged Yukawa fluid (with decay parameter Zq = 2z) as a new 
reference system. Consequently, we are referring to this theory as the FMSA/YO theory. Numerical 
calculations that we performed for the same model discussed by Tang fB^ showed that the proposed 
modifications of the reference system make the first order mean-spherical approximation even more 
accurate. 

It is important that in order to treat the novel reference we are employing the conventional 
FMSA/HS approach. Due to this, the level of simplicity and transparency of the improved 
FMSA/YO theory is kept at the same level as that of the conventional FMSA/HS theory. In 
particular, it is easy to see that the resulting new expressions for main ingredients of the FMSA 
ideology - the Laplace transforms of the radial distribution function, equation ( 17 1, and the Baxter 
factorization function, equation (18 1, of the new reference system, as well as the equation ( [20| for 
the correction term - are only slightly longer than their original counterparts [see equations (|4|, ([5| 
and (13)] and both being composed of the same variables, functions and model parameters. Only 
the parameter of the FMSA/YO theory appears to be the decay parameter Zq, for which in the 
case of Yukawa fluid we obtained Zo = 2z, and which determines an amount of the out-of-core 
attraction that should be attributed to the YO reference system. 
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VflocKOHa/ieHe cepeflMbocctaepuHHe Ha5/iM>KeHH5i nepiuoro nopaflKy 

npOCTMX piflMH 

C. r.nyiiiai^l', A. TpoxMMHyH|il2l_ |. HesGeflaP 

iHCTMTyr 4>i3MKH KOHfleHCOBa H i/ix cwcTeM HAH VKpaTHM, Byji. I. CBeHL4iL4bKoro, 1, 79011 JlbBiB, VKpama 
tpaKy/ibTer xiMiT ra 6ioxiMiV, VHiBepci/iTeT Bpi/iraM 9\Hr, HpoBO, CLUA 
<t>aKy/ibTeT npupoflHUHHx HayK, yHiBepcHTer 9\.E. FlypKiHCie, HecbKa Pecny6jiiKa 

riponoHyeTbcn hobum niflxifl reopiT sSypeHb Ha ochobI cepeflHboc4)epMHHoro Ha6/iM>KeHH!R nepmoro nop^iflKy 
(CCHnn). BiH no/i5irae y Bi/iKopi/icxaHHi TBepfloccfjepHoT (TC) piflUHi/i KDKaBi/i b aKocxi 5a3HCHoT cucreMH, Ha 
ocHOBi nKoT poBB'jiByeTbcn piBHiqHHa OpHLUTeHHa-LlepHiKe b paMKax xeopiT sSypeHb. V poSoxi o6roBopKieTbC5i 
BM6ip onxi/iMa/ibHoro napaMerpa npnTJira/ibhoT fla/ieKOfliT 6a3WCHoT ci/icjeMi/i. Pesy/ibTaTi/i nopiBHKDioTbcn 3i 
BBMHawHOKD Teopieio CCHnn/TC xa MOfle/iKDBaHHiRM MoHTe-Kap/io, fljia KoectDiL4ieHTa CTi/ic/inBocji i cfjaBOBoT 
fliarpaMM piflUHi/i KDKaBH. SanponoHOBana Teopia sBepirae nonepeflHiw piseHb npocroTH Ta npoBopocxi nK i 
BBHHawHe CCHnn/TC, a/ie e 6ijibLU tohhokd. 

K/iiOHOBi c/iosa: pifliAHa K)KaBi/i, cepeflHboc(pepi^HHe Ha6niA>KeHHa, Teopia sBypeHb, MOftenioBaHHa 
MoHTe-Kapno 
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